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The variation in the Arctic sea ice is dom- 
inated by the seasonal cycle with little inter- 
annual correlation. Though the mean sea ice 
area has decreased steadily in the period of 
satellite observations, a dramatic transition 
in the dynamics was initiated with the record 
low September ice area in 2007. The change 
is much more pronounced in the amplitude 
of the seasonal cycle than in the annual mean 
ice area. The shape of the seasonal cycle is 
surprisingly constant for the whole observa- 
tional record despite the general decline. A 
simple explanation, independent of the in- 
creased greenhouse warming, for the shape 
of the seasonal cycle is offered. Thus the dra- 
matic climate change in arctic ice area is seen 
in the amplitude of the cycle and to a lesser 
extend the annual mean and the summer ice 
extend. The reason why the climate change 
is most pronounced in the amplitude is re- 
lated to the rapid reduction in perennial ice 
and thus a thinning of the ice. The analysis 
shows that a tipping point for the arctic ice 
area was crossed in 2007. 

Introduction 

Observations[[I]] show more dramatic decline in 
the ice area as a response to the increased green- 
house gas forcing than predicted by climate model 
simulations [[2 J. Most dramatic is the decrease in 
the September minimum ice extend[[3J [T]]. The 
projections of the Arctic sea ice are important for 
the perspective of future exploitation and possibility 
of commercial shipping through the Arctic Ocean, 
where completely ice free summers have been sug- 
gested as soon as 2030[[3]]. But more importantly, 
the reduction in the sea ice might, through the sum- 
mertime ice albedo feedback, play an essential role in 
the polar amplification of global warming[[5, 6, 7, 8 ]. 
The sea ice extend depends on many factors, such as 
changed atmospheric circulation patterns [[9.], which 
together with changed heat fluxes in the open ocean, 
more clouds and water vapor[[T0]] raise the surface- 
air temperature (SAT) [[5]] and thus increase summer 



melt[[TT]]. Other factors are increased wind driven 
ice export through the Fram Strait [[12]] or changes 
in location and strength of the Beaufort Gyre [[13 J. 
The ice-albedo feedback in the melt season is not 
well understood, it response non-linearly to the 
change in ice ages[ pH [15] ] . Furthermore, the effect 
of the ice albedo also depends on the cloud formation 
over breakup areas [^1$ ]]. The poorly constrained de- 
pendency of the sea ice on all these factors [ [17] ]. 
the large natural variability and the limitations in 
understanding the effect of climate change on pat- 
terns like the North Atlantic Oscillation (NAO) and 
the Northern Annular Mode (NAM) [[9] US]] makes 
projections difficult. Even though there is broad 
consensus that the positive feedback from the ice 
albedo contributes to warming in the Arctic, there 
is much less agreement on it's role in comparison to 
other processes in the polar amplification of global 
warming [[19 ] 120 ] ]. 

In order to better understand how the Arctic sea 
ice has reacted to climate change, an analysis of the 
observed variations in the total ice area is performed. 
The sea ice area is derived from NASAs Satellite 
based Scanning Multichannel Microwave Radiome- 
ter (SSM/I and SMMR)[[3T]] and obtained from 
the University of Illinois' project "The Cryosphere 
Today" [[22]]. The observed record (figure 2, top 
panel), covers the Arctic Ocean and surrounding wa- 
ters in the period 1979-present with daily resolution. 



Decomposing the Annual Cycle 

The resent years, 2007 to present, show a remarkable 
decrease in annual mean (red and green curves in 
figure 1 (a)), with a pronounced summer melting. In 
order to quantify the climate change and change to 
the seasonal cycle the sea ice area can be described 
by 

Xi(t) = mi + Aif(t) + ei{t), 

where i £ (1979, 2011) denotes year, and t £ (1, 365) 
denotes day of year. 
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Figure 1: The seasonal variation of the Arctic sea ice 
area, (a) shows all the years obtained from the satel- 
lite measurements. A year round decline in is seen 
after 2007. The year 1996 showed a notably small 
amplitude in the seasonal cycle, (b) Same data as in 
(a) but normalized by subtracting annual mean and 
dividing by the seasonal amplitude (2012 excluded). 
There is a striking collapse of all years despite the 
pronounced climate change after 2007. (c) The mean 
annual cycle function is obtained as the day-by-day 
average of the 33 normalized curves in (b). 



The function f(t) has zero mean and represents 
the (constant) seasonal cycle. The mean rm and am- 
plitude Ai are constant within a given year i. The 
residual Ei(t) is assumed to be a simple stochastic 
noise. It is not given a priory that this is an ade- 
quate description of the Arctic sea ice area. How- 
ever, by plotting the normalized ice area (xi(t) — 
rrii)/Ai = fit) — €i(t)/Ai we see an almost per- 
fect collapse of the data (figure 1(b)), see the ap- 
pendix for how to estimate m., and Ai. Thus the 
function f(t) can be accurately estimated as the the 
mean of the normalized ice area for all years (fig- 
ure 1(c)). As the difference between the two lower 
panels in figure 1 is small, the residual noise ti(t) 
is small. This is plotted with the ice area in fig- 
ure 2(a) (2012 cannot be included, since the data 
does not cover the full year, however, see supple- 
mentary material: www.gfy.ku.dk/^pditlev/seaice- 
suppl.pdf). Two findings are surprising regarding 
the residual noise: Firstly, despite the dramatic re- 
duction in Arctic sea ice through the record, there 
are no trends in the noise. Secondly, the resid- 
ual noise is a perfect red noise with a correlation 
time of r = 22 days, since the autocorrelation 
is exponential (figure 2(b)) and the compensated 
noise e(t) — e(t) — exp(— l/r)e(t + 1) is perfectly 
uncorrelated and structureless (figure 2(c)). The 
correlation time of 22 days is consistent with the 
time scale of variations in sea ice area governed 
by the natural atmospheric variability. The annual 
mean m, (figure 2(d)), shows, as has been reported 
before[[31[TTJ[531[M]], a downward trend through the 
full record, while the amplitude of the seasonal cycle 
Ai (figurc(2(d)) has a distinct minimum for 1996 and 
a sudden positive jump in 2007, from which it has 
not recovered, indicating a new state dominated by 
one- year ice[[53]]. (See the appendix for a discussion 
on the statistical significance). 



3 
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Figure 2: The decomposition of the Arctic sea ice 
are. (a) shows the satellite measurements since 1979, 
with the small residual noise shown on the same 
scale, (b) The autocorrelation of the residual noise 
is almost perfectly exponential. The green line is 
the curve exp(— t/r), with t = 22 days, (c) The 
scatter plot of the compensated noise (see text for 
explanation) shows that there is no structure, be- 
side the simple exponential autocorrelation in the 
residual noise. (d) shows the steadily declining 
mean ice area rrii and the amplitude of annual cycle, 
which beside the outlier at 1996 shows an increase, 
with a jump in 2007. (e) shows the winter max- 
imum and summer minimum (observe shifted axis 
for comparison), these are obtained from the mean 
and amplitude as mini = m i + Aimin(f(t)) and 
maxi = rrii + Aiinax(f(t)) which occur on Septem- 
ber 8. and March 9., respectively. 



The sea ice area at any particular day of the 
year is easily obtained from the decomposition, 
thus we have the summer minimum mini = mi — 
0.53A, where min(f(t)) — /(September 8) = —0.53 
and the winter maximum maxi = mi + 0.41 A;, 
max(f(t)) = /(March 9) = 0.41. These are shown 
in figure 2(e), both with a downward trend, but not 
as clear a jump in 2007 as is seen in the amplitude 
Ai. Note that the natural variability unrelated to 
climate change, represented by Ci(t) is filtered out of 
the estimates for mini and maxi. 



The Tipping Point 

A change in the linear trend in m, as well as 
in summer minimum can be inferred from lin- 
ear regression analysis, though with low statisti- 
cal significance[|23j. Here we shall implement a 
Bayesian inference change point analysis, which is a 
slightly more advanced method for detecting changes 
in statistics. This analysis, see the appendix for de- 
tails, tests the data for probable changes in statis- 
tics. It is assumed that the data follows gaussian 
distributions with constant mean and variance be- 
tween (unknown) change points. The likelihood of 
the realized data after the previous change point is 
calculated and compared to the likelihood if a change 
of mean and variance has occurred at a new change 
point. This probability is shown in grey scale in fig- 
ure 3(c), where white corresponds to probability 
and black to probability 1. The red curve tracks 
the maximum probability. The contour plot is read 
in the following way: A change point is defined at 
the beginning of the record, thus for a given year, say 
1994, the maximum probability of the duration since 
the last change point is 15 years (time since 1979). 
Probability of duration longer than 15 years is ob- 
viously zero (lower white triangle). A change point 
in the statistics for the amplitude Ai is detected at 
2007. The cumulated distribution function for the 
amplitude prior to the change point is shown in fig- 
ure 3(d), it follows closely a gaussian distribution 
(black curve). Note that the 1996 low amplitude 
is a 3er event, thus in figure 3(c) a slightly higher 
probability of a changing point at 1996 is seen. The 



corresponding analysis for the mean rrii (not shown) 
does not detect any changing points in the statis- 
tics. Thus the amplitude of the seasonal cycle Aj is 
the better fingerprint of Arctic climate change. A 
scatter plot of A, vs. rrii (figure 3(e)) clearly shows 
a change after 2007, which has not recovered (color 
coding is the same as in figure 1(a). Note again 
the 1996 outlier (black)). With the detection of the 
changing point, the mean m, is shown in figure 3(a) 
with tendency for the two periods obtained from lin- 
ear regression (blue lines), in figure 3(b) the ampli- 
tude is shown with the mean values for the two peri- 
ods (blue lines). In both cases the green dashed line 
shows the tendency (for rrii) or mean values (for A4) 
when no change in statistics is assumed. 
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Figure 3: The change in statistics at 2007. (a) 
shows an increase in the linear trend after 2007, 
which showed a deep minimum. Beside the general 
trend (green dashed curve) the change in the slope 
after 2007 (blue lines) is statistically marginal, (b) 
The change in the amplitude in 2007 is much more 
significant. The last five points are much higher than 
the mean for the period (green dashed line) indicat- 
ing a change in mean (blue lines), (c) The result 
of a Bayesian inference of change in statistics of the 
amplitude of the annual cycle. A change point is de- 
tected in 2007 (see text for explanation), (d) shows 
that the annual amplitude prior to 2007, normal- 
ized to unit variance followed a gaussian distribution 
(black curve). Note that 1996 is a 3cr outlier, (e) 
Scatter plot of rrii vs. At where the color coding is 
the same as in figure 1(a). This indicates that 2007 is 
a tipping point to a new statistical state. Note that 
A and rrii are independent within the two popula- 
tions, prior to 2007 (blue and black points) and after 
2007 (red points). 



If all years are considered, there is a significant 
negative correlation between A4 and m, , with a Pear- 
son's r = —0.69. However, if a change is assumed at 
2007, the negative correlation within the two popu- 
lations is much smaller, r = —0.37 (1979-2006) and 
r = -0.26 (2007-2011). In neither case is the cor- 
relation significantly different from zero. This sug- 
gests that the arctic climate has experienced a tip- 
ping point in 2007 to a new state where the mean 
ice area has shifted to a lower level and the ampli- 
tude of the annual cycle has shifted to a higher level. 
The climatic transition is different from the recently 
proposed bifurcation from a mono-stable state to a 
bistable state[ [251l35] ]. 



The Mean Seasonal Cycle 
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Figure 4: The general annual cycle in Arctic sea ice 
area, (a) The cycle shown beginning at the time 
of maximum extent at the dawn of the polar day. 
(b) shows the mean daily insolation at 80N. (c) The 
change in the normalized sea ice area (green) is a 
simple function of the insolation and a melt rate pro- 
portional to the ice area itself. 
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The mean seasonal cycle function f(t) of the ice area, 
is plotted again in figure 4(a) beginning the year with 
the dawn of the polar day. The change in this nor- 
malized ice area can only depend on the wind driven 
transport [[HJ HE]], oceanic flow and SAT in an aver- 
age sense, represented by a constant growth rate Cq. 
Radiative forcing is proportional to the insolation 
I(t) over the arctic ocean and the amount of melting 
ice must be proportional to ice area f(t) itself. Thus 
the simplest possible linear model is: df(t)/dt = 
c - cil(t) - c 2 f(t). The Arctic Oscillation^] and 
other climate fluctuations on time scales longer than 
the correlation time of the noise ej(t) must be re- 
flected in the annual variables rrii and A4. Taking 
I(t) to be the mean daily insolation at 80N (figure 
4(b)), we see that the ice area is at its minimum 
on September 8, right before beginning of the polar 
night, where it grows to its maximum at March 9, 
with the dawn of the polar day. Figure 4(c) shows 
df(t)/dt (in green) and cq — c\I{b) — C2/ (i) (in blue) 
with (c ,ci,c 2 ) = (0.0056, 0.00003, 0.01) day -1 . The 
agreement is surprisingly good. An important find- 
ing is that the accelerations of the ice retreat in the 
polar day and the ice growth at polar night (slope of 
green curve in figure 4(c)) are very similar. This is 
consistent with flux calculations from reanalvsis[|14|] 
that the ice albedo feedback, which is only active 
during polar day is reflected in the the sea ice dy- 
namics through the variables A4 and rrii. 

Summary 

In summary, the Arctic sea ice shows a remarkably 
regular seasonal cycle where only the annual ampli- 
tude and mean change from year to year. This is 
remarkable for two unrelated reasons: Firstly, the 
variations in the total ice area covers regions with 
large variations and different trends [[1 ]. Illustra- 
tive animations are provided by Google Earth[ 28 ]. 
Secondly, different changing factors influencing the 
ice growth and retreat have distinct seasonalities [[17, 
[2"9"] j. The decomposition of the sea ice signal as done 
here, and the constancy of the normalized annual cy- 
cle, indicates that all these factors add up such that 
the amplitude is the strongest indicator of climate 



change. A significant tipping point was crossed in 
2007 where the statistics permanently changed to a 
state of larger seasonal amplitude. 

Appendix: Methods 

Decomposing the sea ice area x(t) 

By averaging equation (1) over the year we estimate 
rrii rs (xi(t)) t ± a/y/h, where a 2 is the variance of 
the noise e(t) and n — 365/22 is the effective num- 
ber of independent points within a year. A poste- 
rior the relative intensity of the noise is calculated 
to ajm = 0.02, thus the uncertainty is negligible. 
Likewise we obtain the amplitude from (xi(t) 2 ) t = 
m 2 +A 2 (f(t) 2 ) t +(T 2 ^Ai = ^(x t (t) 2 ) -m 2 , where 
we have safely neglected the a 2 term, which gives a 
relative error of less than 0.004. The mean cycle 
function is obtained as f(t) = ((xi(t) — rrii)/Ai)i, 
where denotes averaging over the 33 years. The 
uncertainty is of the order a / \{Ai) iV33) = 0.004, 
thus also negligible. As a consistency check of equa- 
tion (1), we obtain (f{t))t/^(f(t) 2 )t = 0.0015 0. 
As to the effect of changing the beginning date for 
the year see the Supplementary material. 

On significance of change in amplitude 

Assuming A^ to be an uncorrelated stochastic vari- 
able, (the autocorrelation is -0.02 for the period 
1979-2006), it is most likely that the arctic climate 
experienced a transition to a a new statistical state 
in 2007. A heuristic argument is based on the obser- 
vation that the five highest values has been observed 
in the last five years of the record. For any distribu- 
tion for Ai, only using the observed independence, 
the probability of a coincidental occurrence of the 
last five consecutive years having the five largest val- 
ues out of 33 years is 1/ (number of ways to pick 5 
years of 33 years)=l/( 3 5 3 ) = 4.2 x 10~ 6 , thus very un- 
likely. Note that this argument for a changing statis- 
tics is not rigorously correct, since any sequence of 
occurrences of, say, the five largest values in 33 draw- 
ings is equally unlikely without indicating change in 
statistics. The fact that the criterion of change is de- 
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rived from the data and not independent flaws the 
argument. Note also that the argument cannot be 
applied to the statistics of the mean m, where the 
value for 2006 is smaller than values for both 2008 
and 2009. 

Change point analysis 

The change point analysis [[30]] is based on Bayesian 
inference, where an unknown number of change 
points in the data series are assumed. In between 
the change points the variables are drawn from a 
gaussian distribution with constant mean and vari- 
ance. The probability of a change point can be cal- 
culated recursively from the conditional probability 
of the next data point given the maximum likelihood 
parameters determining the distribution of the pre- 
vious data. The prior distribution for a change point 
is taken to be uniform. The numerical algorithm is 
taken from ref.[|31]]. 
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